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ABSTRACT 

We present X-ray timing and spectral analyses of simultaneous 150 ks Nuclear Spectroscopic Tele¬ 
scope Array (NuSTAR) and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4i51. We 
disentangle the continuum emission, absorption, and reflection properties of the active galactic nu¬ 
cleus (AGN) by applying inner accretion disk reflection and absorption-dominated models. With a 
time-averaged spectral analysis, we find strong evidence for relativistic reflection from the inner ac¬ 
cretion disk. We find that relativistic emission arises from a highly ionized inner accretion disk with a 
steep emissivity profile, which suggests an intense, compact illuminating source. We find a preliminary, 
near-maximal black hole spin a > 0.9 accounting for statistical and systematic modeling errors. We 
find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry, 
in which the illuminating corona is modeled as a point source. Through a time-resolved spectral anal¬ 
ysis, we find that modest coronal and inner disk reflection flux variation drives the spectral variability 
during the observations. We discuss various physical scenarios for the inner disk reflection model, and 
we find that a compact corona is consistent with the observed features. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: indi¬ 
vidual (NGC 4151) — galaxies: Seyfert — X-rays: galaxies 


1. INTRODUCTION 

The dominant features of the X-ray spectra of Seyfert 
galaxies are a Comptonized continuum, reflection from 
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the accretion disk as well as distant material, and 
line of sight absorption. Uniquely disentangling these 
spectral features has been challenging with most X- 
ray observatories, which primarily have high sensitiv¬ 
ity in the E ^ 0.1 — 10 keV band. High sensitiv¬ 
ity, broadband (0.5 — 79 keV) observations of Seyfert 
galaxies enabled b’ y the Nuclear Spectro scopic Telescope 
Array {NuSTAR) ([Harrison et al.ll201,‘lf l combined with 
observatories with high sensitivity and energy resolu¬ 
tion at lower energies, including Suzaku ([Mitsuda et al.l 
l2007f) . allow the continuum, reflection, and absorption 
features to be definitively separated in Seyfert galaxy 
X-ray spectra. Deconvolving these features can allow 
for fundamental properties of active galactic nucleus 
(AGN) coronae, supermassive black holes (SMBHs), and 
AGN absorption structure s to be determined (e .g., NGG 
1365, IRisaliti et all [Ml. I Walton et all [Ml SWIFT 
,12127.4-1-5654. iMarinucci et al.ll2014a|) . 

The dominant X-ray continuum emission in Seyfert 
galaxies is commonly thought to be produced by Gomp- 
tonization of black hole accretion di sk thermal emission 
in a coronal electron pla sma (e.g. iHaardt et al.l Il99l 
[Reynolds fc Now^ l2003f) . The corona could be pro¬ 
duced by a n accretion disk atm osphere or the base of 
a jet (e.g., iMarkoff et al.l l200ll . Seyfert galaxies can 
produce weak jets ev en though they are radio quiet 
([Ghisellini et al.l 1200411 . For some AGNs, the corona is 
compact with a characteristic distance from the accre¬ 
tion disk D < lOrg (where rg = GM/c^) as inferred from 
X-ray micro-le nsing, reverberation , and eclipse mea¬ 
surements ("e.g.. iG hartas et al.ll200f : iZoghbi et al.ll2?)T2l : 
[Reis fc Milled 1201, 31 : iSanfrutos et al.ll201.3l and references 

therein). Because of the compact sizes inferred for AGN 
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coronae, their detailed geometry has been difficult to de¬ 
termine. 

Coronal emission is reprocessed in the black hole ac¬ 
cretion disk, which primarily produces relativistically 
skewed fluorescent Fe Ka line emission (rest-frame en¬ 
ergy E = 6.4 keV for neutral Fe) and a C ompton scat¬ 
terin g ‘hump’ at energies E > 10 keV fe.g.. lFabian et al.1 
119891 iGeorge fc FabiaDlll99lll . Through measuring the 
relativistic distortion of the inner accretion disk reflec¬ 
tion, important properties of the inner accretion disk and 
black hole can be dete r mined including the black hole 
spin (jBrennemanI 1201,81 iRevnold^ 1201.81 and references 
therein). SMBH spin (which has dimensionless form 
a = cJ/GM^ where J and M are the black hole angular 
momentum and mass, respectiv ely) encodes the accretio n 
and merger history of a SMBH (jBerti &: Volonteiill2008ll . 
Building upon the census of the ^ 22 SMB H spin mea¬ 
surements currently in the literature (e .g., [BrennemaDl 
I201.‘U[Reynoldsl[201.8llWalton et al.ll201.‘ll and references 
therein) is crucial for understanding SMBH eyolution. 
SMBH spin also likely plays a critical role in driving AGN 
outflows and jets t hat deposit in atter and energy in their 
environment le.g.. lFabiaDll20l3l . 

Nuclear emission is absorbed in several r egions com - 
monly associated with ionized outflows fe.g..lElvi iHo^ 
the b road-line region (|Antonnccilll99.8HUrrv fc Padovani 
1995), and/or a molecular torus fe.g. lKrolik fc Begelman 


1988f) . Absorption and relativistic reflection can produce 
similar spectral features in the Fe K band. This has 
led some authors to suggest that the relativistic reflec¬ 
tion features actually arise purely from absorption and 
Compton scattering in material relatively far from the 
black hole where relativistic effects are negligible (e.g., 
iMiller et al.ll200^ . 

The Seyfert 1.5 galax y NGC 4151 (z=0.00332, 
Ide Vaucouleurs et al.l 1199111 . which is sometimes consid¬ 
ered to be the archetypical Seyfert 1 galaxy, is a source 
in which many AG N phenomen a were first characterized 
(for a reyiew, see lUlriciil |2000D . The sou rce has bolo- 
metric luminosity Lboi ~ 5 x lO'^^ergs”^ (|Woo fc Urryl 
[2003), SMBH mass M = 4. 57^^ x 10^ Mq from opti¬ 
cal and UV reverberation (jBentz et al.ll200^ . and cor¬ 
responding Eddington ratio Lboi/LEdd ~ 0.01. Despite 
being one of the brightest and most studied AGN, few 
constraints on the coronal geometry and the black hole 
spin have been found for the source. 

NGG 4151 has a harder X-ray continuum relative to 
the average Seyfert AGN. If approximated with a power- 
law, it has a photon index historically measured to be in 
the range F ^ 1.3—1.9 (e.g. Jives et al.lll976HPerola et al.l 
I986t lY^aaoob fc Warwick 11991 : Zdziarski et al.l Il996t 


Piro et al.Ml999t iBeckmann et ah l2005fl .~ NGC 4151 
was also the first AGN in which a high-en ergy cut- 
off was clearly dete cted (|,Ionrdain et al.l 19921 SIGMA; 
iMaisack et alJ[T9^ CGRO OSSE). From an analysis of 
all INTEGRAL data of the source and data from other 
X-ray obse rvatories during the p eriod 2003 January to 
2009 June, iLubihski et al.l (I2010D found that the coronal 
emission had an approximately constant X-ray spectral 
index and Compton y parameter and inferred that the 
corona had an approximately constant geometry. The 
X-ray continuum conceivably could arise in the base of a 
jet, as NGC 4151 has displayed a non-relativistic jet as 


obser v ed using radio i nterfe r ometry (iWilson fc UlyestadI 
Il982t iMundell et al.l I2003I lUlvestad et al. 12005) and 
chara cterized with optical le.g.. iStorchi-Bergmann et al.l 
l2009fl and X-ray observations le.g.. iWang et al.ll2011al 
Chandra ACIS). 

NGC 4151 has shown inconsistent evidence for re¬ 
flection from the inner accretion disk from X-ray spec¬ 
troscopy. Spectral analyses of observations taken with 
instruments with a E ^ 0.1 — 10 keV bandpass and 
CCD resolution have indicated the pr esence of rel¬ 
ativis tic refl ection in NGC 415 1 (e.g., lYaooob et al.l 
119951 A.S'CA: iNandra et al.ll2007t XMM-NewtonY How- 
ever, other analyse s find no such e vidence for it (e.g. , 
iSchurch et al.l[200l . XMM-Newton; [Patrick et al.ll2oia 
SuzakuESwift). The lack of simultaneous broadband 
[E ^ 1 — 100 keV), high-sensitivity X-ray observa¬ 
tions of NGC 4151 combined with the complex absorp- 
tion (total Nh lO^) ^ cm~^) and soft e mission (e.g., 
lArmentrout et al.1l2007L IWang et al.ll2011bD in the source 
have also made the unambiguous detection of relativistic 
reflection challenging. 

Recent X-ray reverberation studies of eight archival 
XMM-Newton observations have shown that NGC 4151 
has a compact corona surrounding a near-maximally ro¬ 
tating black hole . The discove r y of F e Ka reverberation 
in NGC 4151 bv iZoghbi et al.l ( 20121) indicates the pres¬ 
ence of broad Fe Ka line emission that responds to coro¬ 
nal emission originating from a height above the accre¬ 
tion disk on the order of a few r^. A more re cent Fe Ka 
reverberation analysis bv ICackett et al.l (I2014D shows evi¬ 
dence supporting a maximally spinning (a = 0.998) black 
hole compared to a low spin black hole (a = 0.1) and an 
illuminating source characterized by the lamp post geom¬ 
etry (in which th e corona is modeled as a point source 
on the spin axis; iMartocchia et all I2002L IMiniutti et al.l 
120031) with height h = 7.0jt2)6 errors are given). 

On larger scales, N GC 4151 exhibits partial-covering, 
neutral mat erial fe.g..lHolt et al.lll980D and ionized out¬ 
flows (e.g., IWang et al.l I2011bl and references therein) 
viewed in absorption. The source has displayed ab¬ 
sorption column density variability on time-scales as 
short as At 2 days that indicates that the ab- 
sorbing material is located in the bro ad-line region 
(jPuccetti et al.l[20^lde Rosa et al.ll2007D . The absorp¬ 
tion has also been characterized as having an ionizecP^ 
(^ ^ 500ergcms“^), outflowing [vout 500kms“^) 
compon ent suggestive of an a ccretion disk wind ( Chandra 
HETG: [Kraemer et al.l[^05f) as well as a highly ionized 
^ 10^ — lO^ergcms"^) and blue-shifted {vout ~ 0.1c) 
‘ultra-fast outflow’ (UFO) component that might be as- 
sociated with an ac cretion disk wind or a weak jet 
(jTombesi et al.l[2011l) . 

In this paper, we present timing and spectral analy¬ 
ses of simultaneous 150 ks NuSTAR and Suzaku X-ray 
observations of NGC 4151 in a moderately bright state 
focused on deconvolving the continuum, reflection, and 
absorption features. In a separate paper, we will present 
a reverberation analysis of these observations (Zoghbi 
et ah, in prep.). This paper is structured as follows: 
Sj2] summarizes the observations and data reduction; )J3] 

Ionization is parameterized with the ionization parameter, 
which is defined as the ratio of the ionizing X-ray fiux to the gas 
density: i = LMA, 
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TABLE 1 

Suzaku AND NuSTAR observation information 


Instrument 

Band^ 

( keV) 

Time 

(ks) 

Count Rate 
(ctss“^) 

Counts 

Suzaku XIS-0+3 

2.5 - 9 

150 

6.180 ±0.007 

929099 

Suzaku XIS-1 

2.5 - 7.5 

150 

2.444 ± 0.004 

367408 

Suzaku PIN 

14 - 60 

141 

0.985 ± 0.003 

138881 

NuSTAR FPMAb 

5-79 

141 

5.154 ±0.006 

723597 

NuSTAR FPMB'^ 

5-79 

141 

4.915 ±0.006 

690323 


Note. — Counts and count rates are background subtracted. 

^ Column gives energy bands used in the spectral analysis. 

^ NuSTAR data corresponds to the co-added NuSTAR observa¬ 
tions. 

discusses the timing analysis; gl presents time-averaged 
f ill.II) and time-resolved 1 >14.21) spectral analyses; >jni dis¬ 
cusses systematic errors and the physical implications of 
our results; and >^5] presents our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

Suzaku and NuSTAR observed NGC 4151 from 
November 11-14, 2012. Suzaku observations were taken 
with NGC 4151 placed at the nominal X-ray Imag¬ 
ing Spectrometer (XIS) pointing position. The Suzaku 
observation has observation ID 707024010, and the 
NuSTAR observation is composed of three pointings 
with observation IDs 6000IIII002, 60001111003, and 
60001111005. Suzaku and NuSTAR observation infor¬ 
mation is given in Table [H Source and background spec¬ 
tra provide a source-dominated view of the 2.5 — 79 keV 
emission from NGC 4151 as shown in Figure [TJ We re¬ 
processed all data usi ng HEAsoft ve rsion 6.16, which in¬ 
cludes XSPEC vl2.8.2 (|Arnaudl 1 199611 and xronos v5.22. 

2.1. Suzaku 

We produc e d S uzaku data from the XIS 
(IKovama et all IMll (XIS 0, XIS 1, and XIS 3) 
using the aepipeline script following the Suzaku ABC 
Guid43 with version 22 of the Suzaku software and the 
XIS CALDB release of 2015 January 5, which includes 
updated XIS gain files. We extracted Suzaku XIS data 
products (i.e. spectra, response files, and light curves) 
using xselect and xisresp from circular, 170" radius 
regions centered on the source. For each observation, 
background spectra were extracted from as much of 
the chip as possible. Source and background regions 
excluded all contaminating sources in the field of view 
(primarily the BL Lac object 1E1207.9-1-3945 and the 
LINER NGC 4156) and calibration sources in the 
corners of the chip. 

We rebinned Suzaku XIS spectra and response files 
from 4096 to 2048 channels to expedite spectral fit¬ 
ting without compromising the energy resolution of the 
data. We then co-added data from the XIS front- 
illuminated (XIS-FI) instruments, XIS 0 and XIS 3, us¬ 
ing the addascaspec script to maximize signal-to-noise 
(S/N). We rebinned XIS source spectra into channels 
with a minimum of 25 counts per bin using grppha to 
allow for robust fitting. In the spectral analysis, we 

http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 



Fig. 1.— Source and background spectra for the XIS-FI (shown 
in black), XIS-1 (red), PIN (green), FPMA (blue), and FPMB 
(light blue) instruments. Source count rates are marked with error 
bars, and background count rates are marked with error bars in 
circles. For all instruments, source count rates exceed background 
count rates except for the PIN instrument above E ~ 40 keV. (A 
color version of this figure is available in the online journal.) 

apply a cross-normalization constant to all detectors rela¬ 
tive to the co-added XIS-FI data. XIS-1 fits with a cross¬ 
normalization of 0.96±0.01, which is lower than the nom¬ 
inal XIS-l/XIS-FI cross- normalizatiorPI. 1.019 ± 0.010. 
However, our value is reasonable given that the cross¬ 
normalization depends on the exact spectral shape of the 
source and exactly where the source is on the chip (K. 
Mukai, private communication). 

Because Suzaku and NuSTAR observed NGC 4151 in a 
relatively bright state, we assessed the influence of pileup 
in the XIS data (note that NuSTAR does not suffer from 
pileup) using the pileest script. We found that, for 
XIS 0, 1, and 3, the central ~ 30" has a pileup fraction 
3% fpi ^ 10% assuming a grade migration parameter 
of 0.5. We found that this minor pileup fraction does not 
influence our conclusions in the timing and spectral anal¬ 
yses on the spectral properties of NGC 4151 as described 
in Appendix 

We reprocessed the Suzaku Hard X-ray Detector 
(HXD) PIN instrument (iTakahashi et al.ll2007l) data us¬ 
ing the aepipeline script with the 2011 September 15 
HXD calibration files. We reduced and extracted it fol¬ 
lowing the Suzaku ABC Guid43- We estimated the 
PIN non-X-ray background using the version 2.2 ‘tuned’ 
non-X-ray background (NXB) event file for November 
2012 as well as response and flat-held files for epoch 
11 data. We estimated the PIN cosmic X-ray back¬ 
groun d with XSPEC simulations using the model of lBoldtl 
( 19871). We rebinned PIN data to a S/N of 5 using ISIS 
( Houck fc Denicol^l2000D to facilitate fitting. Three 
percent systematic errors were added to PIN data as 
recommended by the Suzaku team to account for un¬ 
certainties present in the NXB data. We allow the cross¬ 
normalization factor of PIN relative to the XIS-FI instru¬ 
ments to vary in the spectral analysis to improve upon 
the agreement between the PIN and FPMA/FPMB spec¬ 
tra found with the nominal valu43, 1.164 ±0.014. It fits 
to 1.26 ±0.01. 

20 
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http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/watchout.html 
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2.2. NuSTAR 

We reprocessed NuSTAR focal plane module A and 
B (FPMA and FPMB) data using the NuSTAR Data 
Analysis Software (NuSTARDAS) version 1.4.1, using the 
NuSTAR CALDB from 2014 July 1. NuSTARDAS performs 
standard data calibration and data screening, including 
the elimination of data collected when the Earth oc¬ 
culted NGC 4151 and when NuSTAR was in the South 
Atlantic Anomaly. For each of the NuSTAR observa¬ 
tions, we extracted source data products from a 70" 
radius circular region centered on NGC 4151, and we 
extracted background spectra from 70" radius circular 
regions on the same detector chosen to avoid contam¬ 
ination and detector edges. Background regions were 
separated by at least 4' from NGC 4151 and > 7' from 
lE1207.9-b3945 and NGC 4156, which exclude > 98% 
and > 99% of the enclosed coun ts for a point source, re¬ 
spectively (|Harrison et al.l[2013ll . Spectra from the three 
NuSTAR observations were co-added and analyzed sep¬ 
arately for FPMA and FPMB. We rebinned FPMA and 
FPMB spectra into channels with a minimum S/N of 5 
to allow for fitting. We allow the FPMA and FPMB 
cross-normalization constants relative to XlS-Fl to fit 
freely in the spectral analysis, and they fit to 1.04 ± 0.01 
and 1.08 ± 0.01, respectively. 

3. TIMING ANALYSIS 

The Suzaku XIS and NuSTAR light curves and hard¬ 
ness ratios have small but significant variability as shown 
in Figure [2] Notably, a drop in flux coincident with an 
increase in spectral hardness is seen during a At ~ 15ks 
period labelled interval 3. During interval 3, the 3—5 keV 
flux drops by a factor of ^ 1.4 and the 10 — 30 keV flux 
decreases by a factor of 1.2. At the same time, the 
hardness ratio (HR) between the 6 — 10 keV flux and the 
3 — 5 keV flux increases by ^ 5a. Between intervals 4 
and 6 (At ^ 45 ks), the 3 — 5 keV flux rises by a factor of 
~ 1.6, the 10 — 30 keV flux increases by a factor of ~ 1.3, 
and the F{6 — 10 keV)/F(3 — 5 keV) HR decreases by 
~ 4(7. From the beginning of interval 6 to the end of in¬ 
terval 7 (At ^ 40ks), the 3 — 5 keV flux drops by a factor 
of ^ 1.4, the 30 — 79 keV flux decreases by a factor of 
^ 1.3, and the F{6 — 10 keV)/F(3 — 5 keV) HR remains 
approximately constant. 

The root-mean-square fractional variability (RMS 
Fvar) is shown in Figure [31 The overall source variabil¬ 
ity drops dramatically below A ^ 3 keV, where circum- 
nuclear and scattered emissi on dominate the a bsorbed 
emission from the AGN fe.g.. lWang et al.]l2011bll . There 
is also an overall drop in variability with increasing en¬ 
ergy above E 3 keV, which is typical of Seyfert AGN, 
where the more variable model components tend to dom¬ 
inate at softer energies. The sharp, relatively narrow dip 
at 6.4 keV is an indication that the neutral, narrow com¬ 
ponent of the Fe Ka line is not varying much relative to 
the surrounding continuum during the observation. This 
is expected of the neutral, narrow reflection component 
because it likely originates relatively far from the nucleus. 

4. SPEGTRAL ANALYSIS 

We perform time-averaged and time-resolved analy¬ 
ses to characterize the broadband X-ray spectral fea¬ 
tures and variability of NGG 4151. We first undertake 


a time-averaged spectral analysis (i.e., jointly fitting the 
full XIS-FI, XIS-1, PIN, FPMA, and FPMB spectra) to 
understand the components present in the source spec¬ 
trum using the high signal-to-noise (S/N) of the full data. 
We then carry out a time-resolved analysis to character¬ 
ize spectral variability that is averaged out in the full 
observations, which may include variations in the coro¬ 
nal emissio n (which can vary o n hour-long time scales in 
AGN; e.g.. [Markoff et aTll2005ll . the inner accretion disk 
reflection emission ( which can vary on t he order of hours 
in NGC 4151; e.g., iZoghbi et al.ll201^ and absorption 
features ( which can vary on th e order of days in NGC 
4151; e.g.. iPuccetti et aLll2007ll . 

To carry out the spectral analyses, we use XSPEC 
with a standard cosmology of Ho = 70kms~^ M pc~^, 
0.27, and Ha = 0.73 (iKoinatsu et al.ll2011h. We 
use X-ray cross-section values from Verner^et^D ilO^ 
and solar abundance values from iWilms et al.l (j200i 


We include data in the energy ranges indicated in Ta¬ 
ble |TJ 

For the spectral analysis, we apply an inner disk re¬ 
flection model and an absorption-dominated model. In 
brief, the inner disk reflection model inclndes a partial¬ 
covering, neutral absorption component (which we re¬ 
fer to as PCI), a distant reflection component (DRC), a 
cut-off power-law component (PLC), and an inner disk 
reflection component (IDR). The absorption-dominated 
model is identical to the inner disk reflection model ex¬ 
cept that the IDR component is replaced with another 
partial-covering absorber (PC2). For the time-averaged 
analysis, we link all parameters between data sets, and 
we fix redshifts of all model components to the systemic 
redshift. We mode l Milky way absorption with Tbabs 
([Wilms et al.ll2000ll . and we fix the Milky Way absorp¬ 
tion column density to the weig hted-average value fo r 
NGC 4151 from the LAB survey ([Kalberla et al.ll2005fl . 
Nh = 2.3 X 10^°cm“^. We determine 90% confidence 
intervals for one interesting parameter for all parameter 
values presented in the time-averaged and time-resolved 
analyses using the XSPEC error algorithm unless noted 
otherwise. 


4.1. Time-Averaged Analysis 

To develop the time-averaged model, we first charac¬ 
terize the continuum absorption and emission features in 
the 2.5 — 4.0 keV, 7.5 — 10 keV, and 50 — 79 keV bands. 
We model the coronal continuum emission with a cut¬ 
off power-law with photon index F and cut-off fixed at 
Ecut = 1000 keV. Data/model ratios in Figure Sj illus¬ 
trate the presence of a strong absorption cut-off at en¬ 
ergy E ^ 2 — 3 keV indicative of a relatively high line- 
of-sight column density {Nh ^ 10^^ cm“^) as well as a 
prominent, narrow Fe Ka emission line at 6.4 keV and 
Compton hump indicative of reflection. 

We apply a partial-cove ring, neutral absor ption com¬ 
ponent fpartcov*zTbabs: IWilms et al.ll2000il . This ab¬ 
sorber greatly improves the ht (see Table [2]) and sat¬ 
isfactorily accounts for the low-energy continuum ra¬ 
tios as shown in Figure ID We then model the promi¬ 
nent, narrow Fe Ka emission line and Compton hump 
as arising by reflection of cut-off power-law emission 
from Compton-thick, plane-parallel, and relatively neu¬ 
tral (i.e., ^ = lergcms“^) material with the xillver 
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Fig. 2.— Light curves and hardness ratios with 2048 s bins (except for panels (e) and (h), which are shown with 4096 s bins to better show 
the variability of the 30 — 79 keV flux). Note the reduction in the broad-band flux coincident with a general increase in spectral hardness 
seen during a At ~ 15 ks period labelled interval 3. Panels (a) and (b) show the Suzaku XIS-FI 3 — 5 keV and 6 — 10 keV light curves as 
labelled. Panels (c), (d), and (e) show NuSTAR light curves in the 4 — 7 keV, 10 — 30 keV, and 30 — 79 keV bands, respectively. Panels (f), 
(g), and (h) show the hardness ratios F{6 — 10 keV)/-F(3 — 5 keV), -F(10 — 30 'k.eV)/F{A — 7 keV), and _F(30 — 79 keV)/F(10 — 30 keV), 
respectively. Vertical, blue dashed lines indicate time intervals 1-7 analyzed in the time-resolved analysis as numbered. (A color version of 
this figure is available in the online journal.) 
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Fig. 3. — The RMS fractional variability of the data calculated 
with 4000 s bins. XIS-FI is in black, FPMA is red, and FMPB is 
green (XIS-1 is not shown for clarity). Data are included in the 
energy range 0.7 — 10 keV for XIS-FI and 3 — 79 keV for FPMA 
and FPMB. The dashed blue line indicates the rest-frame energy 
of the neutral Fe Ka emission line (6.4 keV). (A color version of 
this figure is available in the online journal.) 


TABLE 2 

AND V AFTER ADDITION OF EACH 
TIME-AVERAGED MODEL COMPONENT 


Additional component 


V 


cut-off power-law 

264344 

4120 

64.2 

partial covering absorber 

15626 

4118 

3.80 

distant reflector^ 

6264 

4112 

1.52 

warm absorber I 

5581 

4110 

1.36 

warm absorber 2 

5513 

4108 

1.34 

inner disk reflector 

4453 

4101 

1.09 


Note. — Each added model component provides 
a significant improvement to the fit. 

^ Values includes a change Ax^/Az/ = —39/ — 4 
from allowing the instrument cross-normalizations 
to fit freely. 

reflection model (|Garcfa et alJl2013ll . We specifically em¬ 
ploy xillver- EcS which allows for a variable cut-off 
energy in the range E = 20 — 1000 keV. We fix the inci¬ 
dent cut-off power-law photon index and cut-off energy 
of this component self-consistently to those of the cut-off 
power-law component, and we additionally fix its incli¬ 
nation to J = 60 deg. As shown in the bottom panels of 
FigureSl the added xillver component significantly re¬ 
duces the ratios towards unity. The remaining residuals 
indicate the presence of blended Fe XXV and Fe XXVI 
Ka {E = 6.67 keV and E = 6.97 keV, respectively) and 
K/3 {E = 7.88 keV and E = 8.25 keV, respectively) ab¬ 
sorption lines as well as broad emission in the Fe K band 
and convex curvature above A ~ 10 keV. 

To account for the unmodeled absorption lines in the 
Fe band, we additionally apply two separate XSTAR ab¬ 
sorption table models, which we refer to as warm ab¬ 
sorbers one (WAl) and two (WA2). WAl is param¬ 
eterized by its column density in the range Nh = 
10^2 _ ]^q 24 ionization parameter ^ in the range 

log(^/ ergcms”^) = 2.5 — 4.5. The absorber has turbu- 

2^ available at http://hea-www.cf a.harvard. edu/~ javier/xillver/ 


lent velocity Vt = 3000 km s~^, which is approp riate for 
a highly ionized absorber (|Tombesi et al.l 120111) . We fix 
the Fe abundance to one. WA2 is parameterized by its 
column density in the range Nh = 10 ^° — 10 ^^cm“^ 
and ^ in the range log(^/ergcms”^) = 0.0 — 4.0. It 
has a turbulent velocity vt = 200kms“^. The Fe abun¬ 
dance is hard-wired at one. We limit ^ to be in the range 
log(^/ergcms”^) = 2.5 — 4.0. For both WAl and WA2, 
^ is sampled with resolution Alog(^/ergcms“^) = 0 . 2 , 
which is important for accurately modeling the absorp¬ 
tion opaci ties and, consequently , determining the black 
hole spin (IRevnolds et al.([2012l) . Having accounted for 
these spectral features, we develop separate inner disk 
reflection and absorption-dominated models to account 
for the remaining broad residuals. 

4.1.1. Inner disk reflection model 

We first characterize the broad residuals in the Fe 
K-band phenomenologically. We begin by applying a 
broad Gaussian. The line fits with equivalent width 
EW = 66 l'l}|eV, a = 0.39 ± 0.04 keV, and line energy 
E = 6.20 ± 0.05 keV. The fit improves by Ax^/Aiz = 
—341/ — 3. The measured line energy and width are 
suggestive of broadened Fe Ka emission. To test for an 
origin from the inner accretion disk, we replace the Gaus¬ 
sian with a model of line emission from the inner accre¬ 
tion disk (relline; iDauser et al.l [ 2 M 0 II with an emis- 
sivity profile e oc r~t with <7 = 3 fixed. The relline 
component fits with line energy E = 6.52 ± 0.01 keV, 
EW = III/I 226 V, inclination i = 6 ± 1 deg, and spin 
a > 0.99. The addition of the broad Fe line is statistically 
motivated by a reduction in by jNv = —354/ —4. 

To account for emission from the inner accretion disk 
self-consistently (i.e., by including the full reflection spec¬ 
trum, including the Compton hump, Fe Ka and E.j5 emis¬ 
sion lines, and the Fe K absorption edge), we include 
an inner disk reflection component (IDR) composed of 
the xillver model con volved with the rel conv relativis¬ 
tic convolution model (jPauser et al.ll20inl l. Specifically, 
we use an optimized version relconvf (Jeremy Sanders, 
private communication), relconv allows the black hole 
spin, a, to fit freely in the range —0.998 < a < 0.998. 
We further allow the inner disk reflection component to 
have a broken power-law emissivity profile that transi¬ 
tions from inner emissivity index <71 to outer emissivity 
index <72 at a radius rbr- relconv additionally assumes 
an inner accretion disk that is geometrically thin, opti¬ 
cally thick, and aligned perpendicular to the spin axis. 

We assume that the inner disk radius (rin) is fixed at 
the innermost stable circular orbit (ISCO) radius and 
that the outer disk radius is fixed at rout = 400 rg. 
This value should be unconstrained assuming 52 ^ 
3 as expected for rap idly spinning black holes (e.g., 
iWilkins fc: Fabianll2011h . We ass ume an isotrop i c emis - 
sion distribution as discussed in iSvoboda et al.l (|2009ll . 
We further assume that the AGN is chemica l ly hom oge¬ 
neous as disc u ssed in, e.g., IRevnolds et al.1 (j2012[ l and 
iWalton et al.l (|2013[ 1. so we fix the Fe abundances of 
the distant reflector and IDR together. We addition¬ 
ally limit the Fe abundance to be in the range 0.5 < 
Ape < 5.0 in order to keep the Fe abundance within the 
range measured for r elatively neutral gas in the source 
(|Yaaoob et al.l 119931 : iSchurch et al.l 120031 : iNandra et al.l 
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Fig. 4.— Top Left: Data/model ratios of a power-law (F = 1.75) absorbed only by gas in the Milky Way that demonstrate the strong 
absorption in the source. Top Right: Data/model ratios after adding a partial covering, neutral absorber and refitting that show the 
spectral features of relatively neutral reflection. Middle Left: Best-fit ratios after additionally including a neutral reflection component 
with free Fe abundance and refitting with F = 1.75 fixed. Middle Right: Same as Middle Left, focused on the ratios in the 4—10 keV band. 
Notice the probable broad Fe Ka line as well as the blended Fe XXV and Fe XXVI absorption lines. Bottom Left: Best-fit ratios after 
additionally including two warm absorber components and refitting with F = 1.75 fixed. Bottom Right: Same as Bottom Left, focused on 
the ratios in the 4—10 keV band. For plotting purposes only, spectral bins have either a significance of 80 (j or a maximum of 80 spectral 
channels for each panel. Also, XIS-FI data is black, XIS-1 is red, PIN is light green, FPMA is dark blue, and FPMB is light blue. (A color 
version of this figure is available in the online journal.) 
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I2007D . We are also motivated to do this to limit the in¬ 
fluence of the degen eracy between the blac k hole spin and 
the Fe abundance (|Revnolds et alJ[2012f) . As shown in 
Table [21 the inclusion of the inner disk reflection compo¬ 
nent is motivated by a reduction in jv by Ay^/Az/ = 
-1060/- 7. 

In summary, our model has the functional form 
(Tbabs)*((partcov*zTbabs)*WA1*WA2*(xillver + 
highecut*zpowerlw + relconv*xillver) ). Best-fit 
parameter values are shown in Table El and the best-fit 
model and data/model ratios are shown in Figure El 

Applying an inner disk reflection model, we find strong 
evidence for relativistic reflection and a preliminary 
SMBH spin constraint a = 0.98 ± 0.01 as shown in Fig¬ 
ure El If we calculate the reflection fraction, R, as the 
ratio of the inner disk reflector flux to the power-law 
component flux in the 20 — 40 keV band, we find that 
R = 1.3 ± 0.2. This reflection fraction is relatively low 
given the near-maximal spin and relatively steep inner 
emissivity index, qi = 10 . 0 /^q' 4 ^, with respect to the max¬ 
imum e xpected reflection fr action in the lamp post ge¬ 
ometry ([Daiiser et al.ll2?)l^ . We discuss this further in 


We find a partial covering fraction that is consistent 
with a scattered fraction. Parameters largely do not 
show degeneracies as illustrated with confidence intervals 
in Figure [7l We do find a degeneracy between the inner 
disk component spin and inclination although we find 
values for these parameters that are consistent with pre¬ 
vious studies as discussed in |EI Foi' ttie cut-off power-law 
component, we find only a lower limit Ecut > 600 keV if 
we allow the cut-off energy to fit freely. For simplicity, 
we hx Ecut = 1000 keV in the best-fit model and find 

p _ 1 yc-f-O.Ol 
i — 1. ZO_g^Q2. 


To illustrate the spectral features that are modeled by 
the inner disk reflection component, we remove this com¬ 
ponent from the best-fit model giving the data/model 
ratios shown in Figure El We then refit this model, 
which gives a significantly worse fit than the best-fit 
model {Ax^/Av = -1-1060/ -I- 7). This increase in 
is driven by residuals in the Fe K-band and the energy 
range A > 10 keV identical to those motivating its in¬ 
clusion. 

We note that non-unity ratios remain at the < 10% 
level in the best-fit data/model ratios at energies E > 
30 keV, as y^/z^ = 1304/1219 (1.07) ignoring data below 
30 keV. We attempt to account for these non-unity ra¬ 
tios with a phenomenological hard power-law component, 
whic h have been observed in r adio-loud AGNs (e.g., 3C 
273; IGrandi fc Palumb^ 1200411 . in order to test whether 
this may b e emission from the dim, extended radio jet in 
the source ([Mundell et al.ll200^IWang et al.ll2011al) . The 
hard power-law component fits with T^ard = 1-231q5q, 
and the coronal power-law P increases to F = I.SSIqq^- 
This component gives a small but significant reduction 
in the fit statistic of Ay^/Az^ = —39/ — 2. The pres¬ 
ence of such a hard power-law component, which could 
perhaps correspond to extended jet emission, would be 
interesting although we caution that the component is 
mainly significant in an energy band where there is more 
uncertainty in the NuSTAR calibration. 

In attempt to constrain the outer extent of the corona, 
we try replacing the relconv component in the best- 


fit model with a relativistic convolution component with 
a twice broken power-la w emissivity profile, kdblurS 
([Wilkins &: Fabiani l2011f) . We allow the inner disk ra¬ 
dius to fit freely because the spin is fixed at a = 0.998. 
The outer index is fixed to = 3 as expected from 
basic geometry at radii in the disk where r ^ h. We 
limit the outer break radius to values r{,r ,2 > 10 fg, and 
we find ri,r ,2 = 10/lop''g- ^nd Vin = 1.5 ± 0.1 rg, 

qi = 10.0+o;4^, rbr,i = 3.5±0.1rg, and ga = This 

model provides only a slightly better fit than the best-fit 
model (Ay^/Az/ = —4/ —1), so a twice-broken power-law 
emissivity profile is not strongly warranted statistically. 

If we app ly _ a lamp post model, relconv_lp 

(jPauser et al.l 120131 1. in place of relconv in the best- 
fit model, we find a lamp post height h = l.d/fg jp rg and 
a reflection fraction consistent with the best-fit model 
(i? = 1.2 ±0.1) although we also find a significantly 
worse fit [Ay^/Au = ±493/ ± 2 relative to the best- 
fit model). We note that the worse fit with relconv_lp 
results from the inner disk emissivity profile being deter¬ 
mined solely with the lamp post height compared to the 
broken power-law profile used in relconv. 

We attempt to add a second lamp post component by 
including a second inner disk component because coro- 
nae could be extended, such as in a scenario where a 
jet produces the coronal emission or a corona has multi¬ 
site activity. We fix all of its parameters to those of the 
first inner disk component except for the normalization 
and lamp post height. We find a similar reflection frac¬ 
tion and an improved fit relative to the single lamp post 
model (Ay^/Azz = —240/ — 2) where the second lamp 
post component has height h = 17 ± 3rg and low flux 
relative to the first lamp post component, which again 
fits with a low height h = l.d/jg'Qp rg. This fit along with 
our best-fit model indicate that a model more compli¬ 
cated than a single point source on the spin axis provides 
a statistically favored description of the reflection spec¬ 
trum. This fit also suggests that the corona may have an 
extended component. 

In attempt to test the extended corona scenario 
self-consistent l y, we apply the relconv±P_ext model 
([Dauser et al.l [201311 . We model the disk-illuminating 
source as an outflow on the spin axis with outflow ve¬ 
locity Vout^ base height hbase, and maximum height htop- 
We fix hbase to the event horizon radius as motivated 
by the dual lamp post model described above. We find 
Vout/c = 0.01/:o;q4p and htop = 1.3/:o;Jprg. This model 
(which gives Ay^/Az/ = ±477/ ± 1 relative to the best- 
fit model) is consistent with the single lamp post model 
presented above within errors. We find a better fit with 
the dual lamp post model than with this extended corona 
model, which may result from describing the corona with 
two physically discrete illuminating components rather 
than a continuous illuminating source. 

We try applying the relxill model in place of the 
IDR component in the best-fit model to more accurately 
account for the emission angle as a function of radius. 
The inner disk reflection models we employ in the best-fit 
model assume a single emission angle for the inner accre¬ 
tion disk emission and also assume an isotropic emission 
distribution. While these are valid assumptions for con¬ 
straining properties of relativistic reflection, the emis- 
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TABLE 3 

Time-averaged spectral analysis best-fit values for the 
INNER disk reflection MODEL 


Component 

Parameter (Units) 

Values 

partial covering 

Nh (cm “2 ) 

(15.4 ±0.6) X 10^2 

(part cov*zTbabs) 

fcov 

0.92 ±0.01 

warm absorber I 

Nh (cm- 2 ) 

2.8/3 3 X 10^2 

(XSTAR grid) 

log(^/ergcms”^ ) 

9 c+O.l 
^•^-O.Op 

warm absorber 2 

Nh (cm -2 ) 

2.4/og X 1022 

(XSTAR grid) 

log(^/ergcms”^ ) 

3.3 ±0.1 

distant reflector 

Xpe 

5.0l°-°^ 

(xillver) 

Kdrc (phcm~ 2 s“l ) 

(1.7 ±0.1) X 10-'^ 

cut-off power-law 

r 

-1 ycr-l-O.Ol 
-*-• ‘ ^- 0.02 


Ecuti keV) 

1000(/) 


EpLC (phcm-^s-l) 

3.7/“;" X 10-2 

inner disk reflection 

a 

0.98 ±0.01 

(relconvf♦xillver) 

<?i 

10.0/°“^ 


<?2 

2.9 ±6.1 


Vbr (l-g) 

3.3 ±0.1 


f'in 

risco (f) 


rout (rg) 

400 (f) 


log(^/ 6 tg cm s~ ^ ) 

3.07 ±0.02 


i (deg) 

IStop 


Kidr (phcm- 2 s-l) 

(1.3 ±0.1) X 10-® 

cross-normalization^ 

XIS-1 

0.96 ±0.01 


PIN 

1.26 ±0.01 


FPMA 

1.04 ±0.01 


FPMB 

1.08 ±0.01 

flux, absorbed'^ 

F 2 _io (ergcm “2 8 “^ ) 

1.50x10-1° 

fit 

X^/i' 

4453/4101(1.09) 


Note. — Statistical errors are quoted to 90% confidence for one 
interesting parameter. The redshifts of all components of the source 
are fixed to the systemic value, 2 : = 0.00332. The Fe abundances 
of the IDR and distant reflector are linked together in the fitting. 
The letter ‘p’ is included in confidence limits that have pegged at a 
limiting value allowed in the model. 

^ Cross-normalizations are relative to XIS-FI. 

^ The model calculated 2 — 10 keV flux is quoted for XIS-FI. 


sion angle of photons leaving the disk is not constant 
and depends upon radius, especially in the innermost 
regions o f the disk as a cons equence of relativistic light¬ 
bending ([Garcia et al.ll20l3l . We find parameter values 
consistent within errors with those for the best-fit model 
including a reflection fraction in the 20 — 40 keV band 
R = 1.3 ± 0.2 and a similar fit to the best-fit model 
(Ay^/Ai/ = —5/ -1-0). Consequently, the assumption 
of a single emission angle for the inner accretion disk 
does not introduce significant systematic error into our 
results. 

4.1.2. Absorption-dominated model 

We investigate an absorption-dominated model 
and develop a model that has the functional form 
(Tbabs) * (WAl * WA2 * (partcov*zTbabs) * 
(partcov*zTbabs) * (xillver + zpowerlw)). We 
find a best-fit model shown in TablelHand FigurelHl This 
model provides a statistically less satisfactory descrip¬ 
tion of the time-averaged data (Ay^/Ai/ = 4-86/ -I- 5) 
relative to the best-fit model. 

The absorption-dominated model also has positive 
residuals at the < 10% level at energies i? > 30 keV, so 
we also attempt adding a hard power-law component to 


it. We find a similar hard power-law component as for 
the inner disk reflection model with Vhard = l-25^g Q 7 
can account for the high energy non-unity ratios. We 
also find a moderate reduction in the fit statistic of 
AxVAi/= -24/-2. 

4.2. Time-Resolved Analysis 

To perform the time-resolved analysis, we produce 
strictly simultaneous time-resolved data products (i.e., 
the data products are time-filtered to only include times 
when the two observatories were simultaneously viewing 
NGC 4151) in seven time-intervals as shown in Figured 
Intervals are chosen to maximize S/N and have approx¬ 
imately constant hardness ratios. We omit making PIN 
time-resolved products due to NuSTAR’s superior sensi¬ 
tivity at £1 > 10 keV. Information for the time-resolved 
data is shown in Table [5j We note that similar time- 
resolved analyses have been performed on joint NuSTAR 
observations with S uzaku or XMM-Newton o f the Seyfert 
galaxies IC 4329A (iBrenneman et al.ll20i^. NGC 1365 

(I Walton et al.l[25Hl . and MCG 6-30-l ^Marinucci et al.l 

I2014bll . 

To first assess the spectral variability in a model inde¬ 
pendent fashion, we analyze difference spectra between 
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Fig. 5. — Top Panel: Model plotted in Ef{E). The components are, from top to bottom, the total model (black), the cut-off power-law 
(dark blue), the inner disk reflection component (red), and the distant reflector component (green). Bottom Panel: The data/model ratios 
for XIS-FI (black), XIS-1 (red), PIN (green), FPMA (dark blue), and FPMB (light blue). The green line indicates a data/model ratio of 
unity. For plotting purposes only, spectral bins have either a significance of 80cr or a maximum of 80 spectral channels. (A color version of 
this figure is available in the online journal.) 
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Fig. 6 .— Ax^ versus the dimensionless spin parameter illus¬ 
trating our statistical spin constraint. The inset shows the range 
a > 0.97. 68%, 90%, and 99% confidence levels for one interesting 
parameter are indicated with dashed lines, which correspond to 
Ax^ values of 1.00, 2.71, and 6.63, respectively. (A color version 
of this figure is available in the online journal.) 


high flux states (intervals 2 and 6) and a low flux state 
(interval 4) as shown in FigureHOl The difference spectra 
peak at E ^ 3 keV and gradually decrease with increas¬ 
ing energy. Additionally, the flux of the narrow, neutral 
Fe Ka line at 6.4 keV varies less than the continuum. 

To fully assess the spectral variability, we perform a 
joint fit of the time-resolved spectra using a strategy 
similar to our time-averaged analysis. We initially set 
the parameter values to their time-averaged model val¬ 
ues (which are shown in Tables [3] and |4]). We require 
many parameters to fit freely to the same value across 
all intervals on physical grounds. We assume the instru¬ 
ment cross-normalizations are constant during the obser¬ 
vations. The distant reflection component, which likely 
arises from material at distances in the range r ^ 7 light- 
days (which corresponds to the approximat e distance to 
the o ptical broad-line region in the source; iBentz et al.l 
[Mil to r ^ 4.2 X10^ light-days (which i s the upper-limit 
on th e size of a dusty torus in the source: iRadomski et al.l 
l2003f) . likely does not vary on the time-scale of the obser¬ 
vations as also suggested by the RMS F^ar of the data. 
The Fe abundance and inclination of the inner accretion 
disk also can be assumed to be constant during the ob¬ 
servations. Similarly, the black hole spin can be safely 
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Fig. 7. — Confidence contours for the inner disk reflection time-averaged model. Confidence contours indicate 99% (black), 90% (red), 
and 68% (green) confidence intervals in two interesting parameters, which correspond to Ax^ values of 9.21, 4.61, and 2.30, respectively. 
Confidence contours are shown for the plane of the PLC F and the IDR normalization (Top Left)^ the plane of a and Ape (Top Right), 
the plane of the PLC F and the IDR ^ (Middle Left), the plane of the IDR a and inclination (Middle Right), the plane of the PLC F and 
the distant reflector normalization (Bottom Left), and the plane of the PLC normalization and the partial covering absorber (Bottom 
Right). (A color version of this figure is available in the online journal.) 
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TABLE 4 

Absorption-dominated model time-averaged best-fit 
VALUES 


C3 

"TD 



Fig. 8. — Data/model ratios are shown after starting with the 
best-fit model and then removing the inner disk reflection com¬ 
ponent without refitting (ratios at top) to illustrate the spectral 
features this component accounts for. Best-fit ratios after refitting 
with r = 1.75 fixed are shown on the bottom. Spectra are rebinned 
strictly for plotting purposes to have bins either with a signal-to- 
noise of SOcr or a maximum of 80 spectral channels. Also, the green 
line indicates a data/model ratio of unity. (A color version of this 
figure is available in the online journal.) 
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Fig. 9. — Top Panel: Absorption-dominated model plotted in 
Ef{E). The broadband components are, from top to bottom, the 
total model (black), the coronal power-law component that is ab¬ 
sorbed by PC2 (red), the power-law component that is not ab¬ 
sorbed by PC2 (blue), and the distant reflector component (green). 
Bottom Panel: The data/model ratios for XIS-FI (black), XIS-1 
(red), PIN (green), FPMA (dark blue), and FPMB (light blue). 
The green line indicates a data/model ratio of unity. For plotting 
purposes only, spectral bins have either a significance of 80 (j or a 
maximum of 80 spectral channels. (A color version of this figure is 
available in the online journal.) 

assumed to be constant, as SMBH angular momentum 
evolves only o n cosmic time-scales from accretion and/or 
mergers fe.g.. iBerti fc Volonteril [20081) . We additionally 
fix the WA2 ionization parameter to its time-averaged 
best-fit value. 

4.2.1. Inner disk reflection model 

We systematically assess which parameters should be 
variable and constant in the time-resolved analysis by 
carrying out fits for all combinations of variable parame¬ 
ters for the following seven parameters: neutral absorber 


Component 

Parameter (Units) 

Value 

PCI 

Nh (cm“2) 

(13.3 ±0.9) X 10^2 


fcov 

0.94 ±0.01 

PC2 

Nh (cm-2) 

56/;4 X 1022 


fcov 

0.42 ± 0.02 

Warm absorber 1 

Nh (cm“2) 

i.etol X io22 


log(^/ergcms“^ ) 

9 C 9 -I-O.O 7 

^•OZ_o 02p 

Warm absorber 2 

Nh (cm~2) 

5 . 9/04 X 1022 


loff(£/ergcms”^ ) 

3.46 ±0.03 

DRC 

.4Fe 

1 c::+0.2 

-*^•^-0.1 


Kdrc (phcm“^s“i) 

( 4.3 ± 0 . 2 ) X 10-^ 


loffff/ergcms”^ ) 

0 (f) 

PLC 

r 

1.92 ±0.01 


Ecut{ keV) 

1000(/) 


KpLC (phcm“2s“l) 

rv -1 9’-|-0.002 

L'-J-J-/ -0.003 

Fit 


4539/4106 (1.11) 


Note. — Statistical errors are quoted to the 90% confidence 
level. Instrument cross-normalizations and model-predicted fluxes 
are consistent with those shown in Table [S] The redshifts of all 
components of the source are fixed to the systemic value, 2 : = 
0.00332. PCI and PC2 indicate partial-covering absorbers 1 and 
2, DRC indicates distant reflection component, and PLC indicates 
power-law component. 


TABLE 5 

Exposure times and background subtracted counts for 
THE Suzaku AND NuSTAR time-resolved data 


Int. 

XIS-FI 

Time Counts 

(s) 

XIS-1 

Counts 

FPMA 

Time Counts 
(s) 

FPMB 

Counts 

1 

2776 

18211 

7245 

2959 

13225 

12946 

2 

19535 

137709 

55138 

18680 

98549 

90109 

3 

6920 

40678 

15810 

6319 

31340 

30462 

4 

24474 

127483 

49964 

23550 

108053 

103254 

5 

7640 

48024 

18871 

7318 

39132 

37702 

6 

2423 

17205 

7043 

2906 

16372 

15494 

7 

6055 

35980 

14411 

7120 

34199 

32830 


Note. — Counts correspond to the 2.5 — 9 keV energy band¬ 
pass for XIS-FI, 2.5 — 7.5 keV for XIS-1, and 5 — 79 keV for 
FPMA and FPMB. The relatively low exposure times for the 
interval 6 data result from relatively poor overlap between the 
NuSTAR and Suzaku observations during this interval. 


Nh, fcov, DRC normalization {Kdrc), PLC P, PLC 
20 - 40 keV flux [Fplc), IDR and IDR 20 - 40 keV 
flux {Flop). We parameterize the normalizations of the 
PLC and IDR in terms of their 20 — 40 keV flux for the 
time-resolved analysis using the cflux model in XSPEC. 
We allow parameters to be variable until we no longer 
find an improvement in the fit based on the F-test. We 
find a best-fit time-resolved fit with four variable param¬ 
eters {Fppc, Fjop, Nh, and P). If we additionally try 
allowing fcov to be variable, it provides a marginally sig¬ 
nificant improvement to the fit according to the F-test 
(Ay^/Aiz = —30/ — 6). We find that this parameter 
is correlated with the power-law normalization, which 
is likely unphysical whether the covering fraction repre¬ 
sents a scattered fraction or a true covering fraction for 
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Fig. 10.— Left: The ratios of the interval 2 (black, top), interval 4 (red, middle), and difference (blue, bottom) spectra to the best-fit 
model to interval 2. Right: The ratios of the interval 6 (black, top), interval 4 (red, middle), and difference (blue, bottom) spectra to the 
best-fit model to interval 6. Notice the narrow dip in both difference spectra at the energy of neutral Fe Ko (6.4 keV as indicated with a 
dashed, light blue line). XIS-FI, XIS-1, FPMA, and FPMB spectra are shown although only NuSTAR data is shown above 9 keV and only 
Suzaku data is below 9 keV for clarity. Spectra have been rebinned for plotting purposes only. The solid green line indicates a data/model 
ratio of unity. (A color version of this figure is available in the online journal.) 


a neutral absorber. Thus, we elect to keep the covering 
fraction constant. 

In summary, the best-fit time-resolved model values are 
shown in Tables m and [T] Time-series of the parameter 
values allowed to vary between intervals (see Figure [TT|) 
show that modest coronal and inner disk reflection flux 
variation drives the spectral variability during the obser¬ 
vations. Note the similarity of the variation in the power- 
law component 20 — 40 keV flux and the 30 — 79 keV flux 
shown in Figure^ The reflection fraction (i.e., the ratio 
of the inner disk reflection and the power-law component 
20 — 40 keV fluxes) is anti-correlated with the power- 
law component 20 — 40 keV flux. This anti-correlation is 
reasonable with respect to some of the predictions of the 
lamp post model for a near-maximal, prograde black hole 
spin. Specifically, for small lamp post heights, the re¬ 
flection fraction can increase when the lamp post moves 
to smaller heights due to a de crease in the continuum 
flux that reaches t he observer (|Miniutti fc Fabi^ 12004 
IDauser et al.ll201^ . 

4.2.2. Absorption-dominated model 

For the absorption-dominated model, we carry out a 
similar systematic analysis as followed for the inner disk 
reflection model. We perform fits for all combinations of 
variable parameters for the following seven parameters: 
the column densities and covering fractions of both neu¬ 
tral absorbers, DRC normalization (Kdrc), PLC F, and 
the PLC 20 — 40 keV flux {Fppc)- We allow parameters 
to be variable until we no longer find an improvement 
in the fit based on the F-test. We find a best-fit time- 
resolved fit with four variable parameters (Fplc? Pj and 
Nh and fcov of PC2) as shown in Figure [121 For the 
time-resolved analysis, we allow these four parameters 
to fit independently in each interval (which are shown in 
Table |5|) and allow all other free parameters (which are 
shown in Table El) to fit freely but linked between each 
interval. 

We find a best-fit absorption-dominated model that 
provides a slightly worse fit of the time-resolved data 


(Ay^/Ajz = -1-29/-1-5) relative to the inner disk reflection 
model. Within the context of the absorption-dominated 
model, variation in the covering fraction of PC2 and the 
PLC flux drives the variability during the observations. 
The PC2 Nh and fcov parameters appear to be corre¬ 
lated with variations in the coronal continuum on the 
short timescales seen in the top panel of Figure |T2l We 
show correlation plots in the bottom panel of Figure [T^l 
to show that fcov and Fppc have closed confidence con¬ 
tours. 

5. DISCUSSION 

Using timing and spectral analyses, we have charac¬ 
terized the constant and variable spectral behavior of 
NCC 4151 as observed simultaneously with NuSTAR and 
Suzaku. Through a spectral analysis, we apply inner ac¬ 
cretion disk reflection and absorption-dominated models. 
We find that the inner disk reflection model provides a 
better fit to the time-averaged and time-resolved spectra. 
Within the context of the inner disk reflection model, we 
find evidence for relativistic reflection from the inner ac¬ 
cretion disk and a preliminary spin measurement for the 
supermassive black hole. 

It is unlikely that a putative neutral absorber (which 
are at radii of r ^ 10"^ — 10^ rg; lAntonuccl I1993L 
lUrrv fc Padovanil 1199511 can partially cover ACN coro- 
nae (which have been found to have si ze D 5, 10 r^ 
for several ACNs including NCC 4151; iReis fc Milled 
l2013f) . Furthermore, there is no physical reason to ex¬ 
pect a neutral absorber to be correlated with variations 
in the coronal continuum on the short timescales seen 
in il4.2.2l because the neutral absorbing gas likely is lo¬ 
cated relatively far from continuum source. A reasonable 
range of distances for the absorbing material might ex¬ 
tend as close as the broad-line region (r = 6.6lJ;g light- 
days; iBentz et al.ll200(t) and as far as the upper-limit on 
the size of a dusty torus in the source (r ^ 4.2 x 10^ 
light-davs: iRadomski et al]l2003[) . In contrast, a corre¬ 
lation between the inner disk reflection and the coronal 
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TABLE 6 

Time-Resolved Best-Fit Values for Parameters Free to Fit within Each Interval 
FOR THE Inner Disk Reflection Model 


Int. 

Nh 

log(F’pic/ergcm 2 s l) 

F 

log(Fjop/ergcm 2 s l) 


(1022 cm-2 ) 




1 


— 9 

-0.09 

1.78 ±0.04 

— 10 14+0-10 
J-U--L^_0.14 

2 


-9.96 ±0.04 

1.77 ±0.02 

— 10 02"*'®'^^ 

1U.UZ_Q 04 

3 


— 10 1 

iU.lZ_Q Q5 

1.73 ±0.03 

-9.88 ±0.05 

4 


-10.22t°-55^ 

1.73 ±0.03 

-9.86 ±0.04 

5 


— 9 QQ+O OS 

y.yo_o.o7 

1 yQ-|-U.UC$ 

-L- <^_0.02 

-9.98t°;°® 

6 


-9-851^8 

1 »7q”I"0.03 
-L- i 'J_o.04 

.-H CN 

d d 
+ 1 
00 
o 
cd 

7 

14.2+5-f 

-W.Obtrol 

1 74+0-02 
-L- <‘±_o.04 

—9 96"^^’^^ 
y.»u_o 06 

Fit 


15655/15589(1.00) 




Note. — 90% confidence intervals are given for the variable parameters for the inner 
disk reflection model, which are the partial covering absorber column density, power-law 
component flux in the 20 — 40 keV band, power-law component photon index, and IDR 
20 — 40 keV flux. 


Interval 


1 2 3 4 5 6 7 



50 100 150 200 250 

Time sinee start of Suzaku observation (ks) 


Fig. 11.— Time-resolved analysis best-fit values for the inner disk reflection model for the variable parameters, which are the partial 
covering absorber column density, power-law component flux in the 20 — 40 keV band, power-law component photon index, and IDR 
20 — 40 keV flux. Vertical error bars indicate 90% confidence intervals, and horizontal error bars indicate the duration of the interval. 
Horizontal lines indicate the mean (solid) and one standard deviation (dashed). (A color version of this figure is available in the online 
journal.) 
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Interval 



50 100 150 200 250 


Time since start of Suzaku observation (ks) 



Fig. 12.— Top: Best-fit values for the time-resolved analysis for the variable parameters for the absorption-dominated model, which are 
the PC2 covering fraction, PC2 column density, power-law component fiux in the 20 — 40 keV band, and power-law component photon 
index. Vertical error bars indicate 90% confidence intervals, and horizontal error bars indicate the duration of the interval. Horizontal 
lines indicate the mean (solid) and one standard deviation (dashed). Bottom: Absorption-dominated model confidence contours for the 
partial covering absorber covering fraction versus the cut-off power-law flux for the time-resolved analysis for each time interval as labelled. 
Confidence contours indicate 99% (black), 90% (red), and 68% (green) confidence intervals for two interesting parameters. 
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TABLE 7 

Time-Resolved Best-Fit Values for Inner Disk 
Reflection Model Parameters Assumed to not Vary 
During the Observations 


Comp. 

Par. (units) 

Values 

PCI 

fcov 

0.93 ±0.01 

WAl 

Nh (cm-2) 

2.9^0° X 10^2 


log(^/ ergcms“^ ) 

2.5(/) 

WA2 

Nh (cm-2 ) 

(2.2 ± 0.9) X 10^2 


log(^/ ergcms“^ ) 

o ^-1-0.2 

DRC 

Ape" 

c n + O-Op 
^•^-0.3 


Kdrc (plicm“2 ) 

(1.6 ±0.2) X 10-4 

IDR 

<?i 

io.o+°;°*> 


<?2 

2 .8t°;? 


rbr (i-g) 

Q q-|-0.2 
'^•‘^-0.1 


a 

0.98 ±0.01 


log(^/ ergcms“^ ) 

3-061^3 


i (deg) 


cross-norm. ^ 

XIS-1 

0.95 ±0.01 


FPMA 

J--UO_o.oi 


FPMB 


Fit 


15655/15589(1.00) 


Note. — 90% confidence intervals are given. WAl and 
WA2 indicate warm absorbers one and two, respectively, PCI 
indicates the partial covering, neutral absorber, DRC indicates 
the distant reflection component, PLC indicates the power- 
law component, and IDR indicates the inner disk reflection 
component. 

The Fe abundances of the IDR and DRC components are 
linked together in the fitting. 

^ Cross-normalizations are relative to XIS-FI. 

continuum is physically possible based on light crossing 
time arguments and has been observed. T he distance be¬ 
tween the corona and the disk inferred bv iCackett et al.l 
(|2014fl and the size scale of the inner disk, which are both 
< 10 rg, correspond to a light-travel time < 2 ks assuming 
the black hole mass M = 4.57lo'47 x 10^ M© (|Bentz et al.l 
1200611 . The reverberation analysis of iZoghbi et al.l (I2012D 
shows that the inner disk reflection flux can lag behind 
the coronal flux on time-scales ^ 2ks. 

Thus, on statistical and physical grounds, we discuss 
only the inner disk reflection model. Here, we discuss 
assumptions that may introduce systematic errors in the 
preliminary black hole spin value and physical scenar¬ 
ios for the disk-corona system. For calculations that use 
the black hole mass, we use M = 4.57^0 47 x 10^ Mq 
(|Bentz et al.ll2006l l measured from optical and UV rever¬ 
beration. 

5.1. Systematic Uncertainties in the Black Hole Spin 

Assumptions in our model introduce systematic uncer¬ 
tainty into the derived confidence intervals in the dimen¬ 
sionless black hole spin. A primary assumption of mod¬ 
eling the inner disk reflection with the relconv kernel 
is that the accretion disk is radiatively efficient, opti¬ 
cally thick, and geometrically thin at radii greater than 
or equal to the ISCO radius. These assumptions gener¬ 
ally should be applicable to NGC 4151 considering the 
modest accretion rate implied by its bolometric luminos¬ 
ity Lboi ~ O.OlLEdd- The assumption of a geometrically 
thin disk breaks down where the disk thickness, t, is on 


the same order as the radius in the disk, r, or t/r ^ 1. 
For Seyfert galaxies with Eddington ratios in the range 
Lboi/LEdd ~ 0.01 — 0.3, we e xpect disk thicknesses in the 
range t ^ (0.1 — 1.0) I'g for a IShakura &: SunvaevI (Il973f) 
disk. Thus, this assumption may be invalid for Vm ^ 2 rg, 
which corresponds to spins a > 0.94. Thus, considering 
the finite thickness of the disk close to the black hole, we 
estimate that this assumption relaxes our spin constraint 
to a > 0.94. 

We assume that the inner accretion disk radius is the 
ISCO radius and that a negligible level of emission orig¬ 
inates within that radius. The latter assumption is sup¬ 
ported by the near-maximal spin and relatively low Ed¬ 
dington ratio of the source, which may indicate that some 
of the accretion flow within the ISCO radius is optically 
thin (|Revnolds fc Fabianll200^ . Additionally, when we 
allow the inner edge of the accretion disk to fit freely, 
it fits to a value consistent with the ISCO radius al¬ 
though we restrict it t o fit to values Tin > rrscn . Based 
on MHD simulations (iRevnolds fc Fabianll200^ . if sig¬ 
nificant emission comes from within the ISCO radius, 
the fit will be biased to higher spins systematically with 
error < 2% for rapidly spinning, prograde black holes 
(a > 0.9). Thus, we estimate this assumption adds ~ 2% 
systematic uncertainty to our spin value. Including the 
systematic uncertainties introduced by our assumptions 
about the thickness and extent of the inner disk, we find 
a constraint a > 0.9. 

Black hole spin has some degeneracy with the iron 
abundance of the inner disk, as both properties affect 
the s trength of the broad Fe line (e.g., IRevnolds et al.l 
l2012f) . We find that NGC 4151 has a super-solar Fe 
abundance. There is no apparent degeneracy between 
our derived Fe abundance and spin values as shown in 
Figure [71 Fe abundance is constant on long time-scales. 
There is evidence for super-solar Fe abundance in mate¬ 
rial distan t from the SMBH. F rom several Ginga mea¬ 
surements, lYaaoob et al.l (I1993D report an average value 
Ape ~ 2.5 from modeling an Fe K absorption edge. From 
analyzing the strength of a neutral Fe abso r ption edge in 
XMM-Newton observations. ISchurch et ^ (120031) report 
an Fe abundance Ape ^3.6, and they report Ape ~ 2.0 if 
they associate the Fe abundance of the neutral absorber 
with that of a Compton reflection compon ent. _ 

For the Fe abundance of the inner disk, iNandra et al.l 
(|2007fl measure an Fe abundance in the range Ape = 
0.3 — 0.5 from analyzing three XMM-Newton observa¬ 
tions. They assume the distant reflection component and 
inner disk reflection component (which are modeled with 
pexmon and kdblur2*pexmoii, respectively) are chemi¬ 
cally h omogeneous in the s ame way as in the present 
work. INandra et al.l (120071) caution that their XMM- 
Newton values require high sensitivity, broad-band ob¬ 
servations that can simultaneously constrain the Fe Ka 
line and Compton hump reflection features. Our analysis 
simultaneously constrains both features, and our derived 
Fe abundance (Ape = 5.0lg;°^) is moderately high in 
comparison with measurements of super-solar Fe abun¬ 
dance in material distant from the SMBH. This may in¬ 
troduce some bias in our spin constraint as a result of 
the po sitive correlation bet ween spin and the Fe abun¬ 
dance ([Reynolds et al.ll2012n . However, we conclude that 
our derived Fe abundance negligibly biases our spin con- 
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TABLE 8 

Time-Resolved Best-Fit Values eor Parameters Free to Fit within Each 
Interval eor the Absorption-Dominated Model 


Int. 

fcov 

Nfi 

logfEpuc/ergcm l) 

r 



(10^^ cm“2 ) 



1 

0.30 ±0.08 

55± 16 

-9.86t°;“ 

1 94+0-07 

l-y^_0.06 

2 

0.35 ±0.03 

64± 10 

- 9 - 78 + 0.22 

1.95 ±0.03 

3 

u.4y_o.o4 

6 + 7 

-9.78l2:°i 

1.93 ±0.04 

4 

Q Cl +0.02 

^•DJ^_0.03 

57l® 

-9.8+2:°i 

1 Ql+O-Oi^ 

i-yi_o.o3 

5 

U-+^_0.05 

40 ±5 

—Q 74+0.01 

y+^_o.o2 

1.83 ±0.03 

6 

0 32+0.07 

U.OZ_Q 08 

47± 13 

-9.73t2:°i 

1 Qc:+y-0t> 

J--oo_o.o6 

7 

0.41 ±0.05 

61+ 

-9-8112:2? 

1.91 ±0.04 

Fit 

to 

15684/15594(1.01) 



Note. — 90% confidence intervals are given for the variable parameters for the 
absorption-dominated model, which are the PC2 covering fraction, PC2 column den¬ 
sity, power-law component flux in the 20 — 40 keV band, and power-law component 
photon index. 


TABLE 9 

Time-Resolved Best-Fit Values eor Parameters 
Assumed to not Vary During the Observations for the 
Absorption-Dominated Model 


Comp. 

Par. (units) 

Values 

PCI 

Nh (cm“^ ) 

13.5l?'2 X 10^^ 


fcov 

0.94 ±0.01 

WAl 

Nh (cm-2) 

(1.6 ± 0.5) X 10^2 


log(£/ergcms“^ ) 

2-5(/) 

WA2 

Nh ( cm-2 ) 

5.912 ? X 1022 


los(£/ergcms”^ ) 

3.5 ±0.1 

DRC 

Ape" 

1 '7+0.2 
-*-•'-0.3 


Kdrc (phcm-2 s-l ) 

(4.1 ±0.3) X 10-4 

cross-norm. ^ 

XIS-1 

0.95 ±0.01 


FPMA 

1.03 ±0.01 


FPMB 

1.06 ±0.01 

Fit 

X^/u 

15684/15594(1.01) 


Note. — 90% confidence intervals are given. WAl and 
WA2 indicate warm absorbers one and two, respectively, PCI 
indicates the partial covering, neutral absorber, DRC indi¬ 
cates the distant reflection component, and PLC indicates the 
power-law component. 

The Fe abundances of the IDR and DRC components are 
linked together in the fitting. 

Cross-normalizations are relative to XIS-FI. 

straint of a > 0.9. 

Black hole spin can show a degeneracy with the disk 
inclination if the blue wing of the broad Fe line profile is 
not significantly detected, as both parameters affect the 
width of the broad Fe line profile. A degeneracy is seen 
between the disk inclination and spin shown in Figure [71 
The inclination of the inner disk should be constant on 
long time-scales, so we compare our determined incli¬ 
nation angle with various inclination angles previously 
reported for the source. 

The inner disk inclination angle we find is consis¬ 
tent with most values determined for the inner disk. 
From the ana l ysis o f three XMM-Newton observations, 
iNandra et al.l (l2007fl find an inclination of the inner disk 
of i = 17^17 deg, * = 21^21 deg, and i = SSlgdeg, 
respectively, where values are given for the best-fit 


kdblur 2 *pexmon model prese nted for each o bserv ation. 
From an ASCA observation, lYaaoob et al.l (|1995f) find 
i = Oito^ deg using the diskline model ([Fabian et al.l 
1198911 . iCackett et al.l (1201411 report i < 30 deg for the 
inner disk from Fe Ka reverberation measurements. 

Comparing the inner disk inclination angle to that of 
more distant structures, our determined inner disk incli¬ 
nation angle is not consistent with the inclination of the 
narrow-line region, 45 ± 5 deg, as determined from Hub¬ 
ble Space Telescope observations (iDas et al.ll2005ll . This 
suggests that the inner disk and narrow-line region are 
misaligned. The inclination we find is consistent with 
the galactic disk inc lination^ i 2 1 deg, as inferred from 
optical photometry (|Simkinlll975[l and Very Large Ar ray 
observations of neutral hydrogen ([Pedlar et al.lll99^ . 

It is possible that some of the broad Fe K line pro¬ 
file may arise from an alternate mechanism than reflec¬ 
tion from the inner accretion disk. For example, scatter¬ 
ing in a Compton thick win d near the disk can produce 
broadened Fe line profiles (|Sim et al.l 120081 1 . However, 
the strength of this feature significantly depends upon 
the outflow mass loss rate. A low mass outflow rate is 
expected for the low mass accretion rate of the source im¬ 
plied by its Eddington ratio Aboi/Afidd ~ 0.01. Thus, we 
estimate that any contribution of scattering in a Comp¬ 
ton thick wind near the accretion disk to the broad Fe 
line is relatively weak compared to the broad line profile 
from inner accretion disk reflection. 

5.2. Physical Scenarios for the Corona and Inner 
Accretion Disk 

Our derived high inner disk emissivity index {qi = 
10 . 0 ^ 04 ^) indicates that the corona that illuminates the 
inner disk is compact. The inner emissivity index and 
photon index, F = 1 . 75 to;o 2 ) together indicate a lamp 
post height h < 10 rg assuming a lamp post geometry, 
as the maximum inn er index for F 1 . 8 , h = 10 rg, and 
a = 0.99 is ~ 3 ([Danser et al.l 1201^ . A lamp post 
height h < 10 rg is expected in order to produce strong 
reflection featur es from which black hole spin can b e 
well constrained ([Danser et al.ll201^ iFabian et al.ll201^ . 
Furthermore, the high ionization parameter we find for 
the inner disk atmosphere is expected for a close illu- 
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minating source. We note that a high emissivity in¬ 
dex {qi > 3) is only physically consistent with a near- 
maximal black hole spin. 

We calculat e theoretical emissivi ty profiles for the lamp 
post model of iDauser et al.l (1201311 for our best-fit F and 
a values. For this calculation, we assume a minimum 
height in order to give a rough estimate of the highest 
emissivity index expected. We estimate a value of ~ 7 
if we average the emissivity from rjsco to the emissivity 
break radius that we find for the best-fit model (r^^.). So, 
we find an inner emissivity t hat is steeper t han ex pected 
for the lamp post model of iDanser et al.l (120131 1. This 
partially illustrates why the broken power-law emissivity 
profile for the inner disk provides a better fit to the data 
than the lamp post emissivity profile. However, this esti¬ 
mate does not necessarily indicate that our best-fit model 
is unphysical. 

We find a relatively low reflection fraction compared 
to the maximum expected reflection fraction, which is 
i? 5 — 10 for the lamp post model descr ibed in M.l.ll 
with a = 0.98 ± 0.01 and h = l.d^o^oprg (jPauser et al.l 
1201411 . This value is a factor ^4 — 8 higher than what 
we find. There are several possible physical scenarios for 
the corona and inner accretion disk that may cause this 
apparent discrepancy. 

If the corona results from a mildly relativistic outflow 
rather than a static region, then the reflection fraction is 
expec ted to be reduced because of relativisti c aberration 
(e.g., iBeloborodovlllOO^ iMalzac et al.ll200lll . A photon 
index T ~ 1.75 indicates v/c ^ 0.4 near the base of an 
AGN jet and a reduction in the strength of reflection, 
H/27r, by a factor of ^ 10 using the disk inclination i = 
18 deg (|Beloborodovlll99^ . This can help explain the 
apparently low reflection fraction we find. We note that 
VLBl observations have shown that the radio jet in NGC 
4151 has had a component with velocit y v < 0.050c at 
a dis tance r ^ 0.16 pc from the nucleus (|Ulvestad et al.l 
1200511 . Thus, if the corona corresponds to the base of a 
mildly relativistic jet, then a mechanism is required to 
decelerate the jet to non-relativistic speeds on ~ 0.1 pc 
scales. Also, an outflow t ends to flatten the emissivity 
profile (|Dauser et al.ll201^ . Keeping these facts in mind, 
an outflowing corona may help explain the low reflection 
fraction we find. 

Truncation of the inner edge of the accretion disk at 
radius r > rrsn n can similarly reduc e the expected re¬ 
flection fraction. iLubihski et al.l (j2010t l conclude that the 
relatively weak strength of reflection from the disk they 
measure (r2/27r ~ 0.3) from all INTEGRAL observations 
of NGG 4151 from 2003 January to 2009 June can be 
explained by an inner hot accretion flow surrounded by 
a cold disk truncated at ~ 15 rg or, alternatively, a 
mildly relativistic coronal outflow. When we allow the 
inner radius of the accretion disk to fit freely using the 
kdblurS model as described previously, it fits to a value 
consistent with the IS GO radius of a black hole with the 
near-maximal, prograde spin that we find for the best fit 
model. Thus, we do not find evidence for a truncated 
disk from spectral fitting of our observations. 

Dilution of the reflection features may also reduce 
the measured reflection fraction. For a close disk- 
illuminating source suggested by the steep emissivity pro¬ 
file, it is plausible that the inner disk sees an intense ir¬ 


radiating flux and has a tenuous, hot atmosphere. This 
may have a density profile roughly similar to that of an 
accretion disk in hydrostatic equilibrium, for which an 
intens e irradiating flux can lead to disk thermal insta - 
bility (jNavakshin fc Do^ 120011 : iBallantvne et al.ll2001fl . 
This instability can produce a highly ionized ‘skin’, 
which could both further broaden the disk reflection 
features and reduce th e observed reflection fraction. 
IBallantvne et al.l (l2001f) showed that simulated RXTE 
and ASGA data of hy drostatic reflection models (e.g., 
IBallantv ne et al.l 120011) fit with c onstant density mod- 
els (|Magdziarz fc ZdziarsO 119951: iRoss fc FabianI [r993D 
tend to fit with systematically low reflection fraction val¬ 
ues because of dilution. Although NuSTAR provides 
superior sensitivity in the 20 — 40 keV energy band, 
where the reflection spectrum is almost completely in¬ 
dependent of ionization parameter, our fit is biased to¬ 
wards best matching the spectral shape at lower energies 
{E < 20 keV), where the reflection spectrum is sensitive 
to the ionization parameter. 

With the xillver models, we assume that the disk is 
illuminated at 45 deg. However, for a compact coronal 
geometry, the disk illumination may be incident on the 
disk at closer to gr azing angles, especi ally within a few rg 
of the black hole (iDauser et al.l[2M^ . In the innermost 
regions of the disk, this would enhance the influence of 
dilution and thus further reduce the reflection fraction. 

In summary, both dilution from a highly ionized disk 
skin and a mildly relativistic outflow may contribute to 
the reduced reflection fraction that we find relative to 
predictions for the lamp post geometry. Because an out¬ 
flow tends to flatten the emissivity profile, the forma¬ 
tion of a highly ionized disk ‘skin’ may be the dominant 
cause of the reduced reflection fraction. We do not find 
evidence for truncation of the inner accretion disk. 

5.3. Coronal Properties 

We model the coronal emission with a cut-off power- 
law with Ecut = 1000 keV fixed and find F = 1.75lQg2- 
This value is consistent within 3a with the average 
Seyfert photon index, F = 1.84 ± 0.03, measured for 
105 Se yfert ga l axies with redshift z < 0.1 with Bep- 
posax (iDadinal [200^1 . It is relatively hard compared 
to the average photon index, F = 1.93 ± 0.01, for 144 
Seyfer t galaxies in the secon d INTEGRAL AGN cata¬ 
logue (jBeckmann et al.l 12009( 1. The value we find agrees 
within ItT with the correlation between the photon in¬ 
dex and Eddington ratio for a sample of 6 9 radio-quiet 
AGNs out to z ~ 2 (jBrightman et al]l2013f) . which gives 
F = 1.63 ± 0.16 for Lboi/LEdd = 0.01. We note that the 
significant dispersion in this correlation as well as uncer¬ 
tainty in estimating the Eddington ratio should be noted 
when applied to a single source. 

The photon index we find is consistent with the range 
in pho ton index (F ~ 1.7 — 1.86) that ILubihski et al.l 
(|2010f) measure for the source in a large range of X- 
ray flux states as observed with INTEGRAL during the 
period 2003 January to 2009 June. The energy cut-off 
we find is corr o borat ed by the energy cut-offs found by 
ILubihski et al.l (I2010D for the source in high and medium 
flux states {Ecut = 264126 keV and Ecut > 1025 keV, 
respectively), which serve as upper and lower bounds, re¬ 
spectively, to the 20 — 100 keV flux of the source in our 
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data. 

6. CONCLUSIONS 

We present X-ray timing and spectral analyses of con¬ 
temporaneous, 150 ks NuSTAR and Suzaku X-ray ob¬ 
servations of the Seyfert 1.5 galaxy NGC 4151. We test 
separate inner disk reflection and absorption-dominated 
models in order to determine the properties of the inner¬ 
most regions of the AGN. Our primary conclusions are 
enumerated below: 

1. We find that the inner disk reflection model provides 
a better statistical and physical description of the data 
than the absorption-dominated model for both the time- 
averaged and time-resolved spectral analyses. 

2. Within the context of this inner disk reflection 
model, we find that relativistic emission originates in a 
highly ionized inner accretion disk with a steep inner 
emissivity profile, which suggests a bright, compact in¬ 
ner disk-illuminating source. We find a relatively mod¬ 
erate reflection fraction with respect to predictions for 
the lamp post geometry. We find a preliminary, near- 
maximal black hole spin a > 0.9 accounting for statistical 
and systematic modeling errors. 

3. Through a time-resolved spectral analysis, we find 
that modest coronal and inner disk reflection variation 
drive the spectral variability during the observations. 

4. We discuss physical scenarios with respect to our 
derived inner accretion disk and coronal properties. We 
find that a compact coronal geometry can reproduce the 
observed features. 

This work further displays the ability of NuSTAR in 


conjunction with X-ray observatories such as Suzaku., 
XMM-Newton, Chandra, and Swift to robustly constrain 
coronal and accretion disk properties in Seyfert galaxies. 
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APPENDIX 

INFLUENCE OF THE MINOR PILEUP IN SUZAKU XIS DATA 

To investigate the influence of the mild pileup in the XIS data on our results, we extract XIS data from an annular 
region with inner radius 30" and outer radius 170". We fit our best-fit model to this data, and we find consistent 
parameter values within 90% confidence in one interesting parameter except for instrument cross-normalizations relative 
to XIS-FI. We find cross-normalizations for XIS-1 of 0.93±0.01, PIN of 1.21±0.01, FPMA of l.OOiO.Ol, and FPMB of 
1.03 ±0.01. Thus, we find that the minor pileup fraction does not influence our conclusions on the spectral properties 
of NGC 4151. 
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